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We study two geometrical factors needed for the correct construction of statistical en- 
sembles of surfaces. Such ensembles appear in the study of fluid bilayer membranes, though 
our results are more generally applicable. The naive functional measure over height fluctu- 
ations must be corrected by these factors in order to give correct, self-consistent formulas 
for the free energy and correlation functions of the height. While one of these corrections — 



P3 ' the Faddeev-Popov determinant — has been studied extensively, our derivation proceeds 

O . 

CJ ■ from very simple geometrical ideas, which we hope removes some of its mystery. The other 

factor is similar to the Liouville correction in string theory. Since our formulas differ from 
those of previous authors, we include some explicit calculations of the effective frame ten- 
sion and two-point function to show that our version indeed secures coordinate-invariance 
and consistency to lowest nontrivial order in a temperature expansion. 
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1. Introduction and Summary 

Amphiphilic molecules in water tend to aggregate into thin flexible bilayers, which in 
turn form a wide variety of nearly two-dimensional structures with characteristic size of 
order microns [1]. The study of these structures may become significant for biology: they 
are certainly important technologically as the underlying elements of microemulsions and 
other complex fiuids. 

To understand the static and dynamic properties of bilayer membranes, and the tran- 
sitions between different morphologies, we must first understand the role of thermal fiuc- 
tuations in two-dimensional surfaces, a subtle problem in statistical mechanics. Thermal 
fiuctuations are important because regardless how stiff a membrane may be, very long- 
wavelength undulations in its shape are always allowed; on long enough length scales any 
membrane will appear flexible and will undulate significantly, a phenomenon first described 
in the cell walls of red blood cells in the 19th century [2].^ 

Since the length scales we wish to study are much larger than the size of the constituent 
molecules, we expect that a membrane will be characterized by just a few effective param- 
eters, summarizing the effects of the complicated molecular forces. Indeed, the famous 
Canham-Helfrich model[4] describes the equilibrium statistical mechanics of membranes in 
terms of just two parameters: a chemical potential for the addition of surfactant molecules, 
which we will call /Ug, and a bending energy coefficient, the stiffness Kq. (Another param- 
eter, the gaussian stiffness Kg, only enters when we consider topology change.) The bare 
coefficients /Uq, Kq are not directly observable, but from them we can derive some more rel- 
evant phenomenological parameters. We will mainly study the "frame tension" r, which is 
the free energy per unit area of the fiuctuating membrane, and the two-point correlation 
{h{u)h{0)) of the height h{u) of the membrane from its average plane, which we fix by 
fixing the edges of the membrane to lie on a square frame. The leading behavior of this 
two-point function at small wavenumber defines another phenomenological parameter, the 
"g^ coefficient" r via 

where A^ is the area of the frame, T is the temperature (we set Boltzmann's constant k^ 
equal to one), and g is a wavenumber. Another reason to introduce r is that in dynamics 



Strictly speaking we are here talking about "fluid membranes," those with no internal struc- 
ture or order within the surface. In this paper we discuss only fluid membranes, but our measure 
factors can be used in more complicated situations too [3]. 
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problems we can still calculate correlation functions like (1.1), while there is no analog of 
the free energy from which to compute r. Two of us study the dynamical problem in [5]. 

In this paper we will explore the relation between the various coefficients /Uq, Kq, r, 
and r. More generally, on long length scales the precise size of the constituent molecules 
is immaterial: we can get the same answers using molecules of size o with coefficients Hq, 
Kq or with rescaled molecules of size 6~^a, 6 < 1, and effective coefficients /Ugj^ (6), K,^g{b). 
The scale dependence of Kgg is well known [6]-[ll], but various answers for (Ugg have been 
given (for example, see [8][9][11]-[15]). In part the differences reflect convention, but there 
is a key physics point which we will address: to get correct answers we must define our 
functional measure properly. Two of us have already discussed such matters in [3], but 
here we will add a number of points, as follows. 

We first introduce the effective action r[h] of the fluctuating membrane and recall the 
argument of [3] that its form is constrained by the coordinate invariance of the underlying 
theory. In particular we show in sect. 2 that the coefficients r and r defined above must be 
equal. We then set up a naive, uncorrected, calculational scheme which yields (consistent 
with [15]) 
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where A = 2n/a and a is the linear size of the constituent molecules. Following Morse and 
Milner we have introduced the thermal de Broglie wavelength 



A = h/V2nmT , (1.3) 

where m is the mass of a surfactant molecule [16]. The same calculational scheme however 
yields a very different formula for r. In fact (1.2) is our final result, but clearly we need to 
work a bit harder to understand r. 

The first correction factor, the "Faddeev-Popov" factor, is well known. In [3] we gave 
a formula for this factor which differs from the one given by David [14]; here in sect. 3 
we give a simple geometrical motivation followed by a derivation of our formula which is 
simpler than the one in [3]. While unimportant in the calculation of r, this factor does 
modify the naive calculation of r (and higher correlation functions too) , yielding an rpp 
which still disagrees with r. 

Next in sect. 4 we argue that another correction factor is needed. Physically our 
membrane consists of constituent molecules each occupying a fixed area a^ in physical 
3-space. As the surface bows outward away from its flat equilibrium configuration its area 



increases and additional molecules must be brought in from solution, effectively increasing 
the number of degrees of freedom in the membrane problem. A statistical measure with 
variable number of degrees of freedom, depending on the configuration itself, is rather 
complicated. We would prefer a measure with fixed number of degrees of freedom, as indeed 
we used implicitly in the naive derivation of rpp. Since changing the number of degrees of 
freedom, or equivalently the ultraviolet cutoff, can be compensated by renormalizing the 
bare energy, we expect that a counterterm correction will need to be added to the naive 
derivation. 

We fix this counterterm by requiring consistency between two different Monge-gauge 
calculations of the free energy (in the Appendix we give a different strategy). The required 
counterterm turns out to make a contribution precisely changing the q^ coefficient from 
Tpp to r = r. 

The introduction of our second correction factor may seem like ad hoc wish- fulfillment, 
but in the Appendix we give some detailed calculations showing how it and the Faddeev- 
Popov factor together are crucial to obtain consistent answers; in particular we show that 
the specific coefficient of the counterterm implied by consistency is the right one to get 
r = T. Furthermore, an analogous correction is well known in the string theory literature, 
where it is called the "Liouville factor" [17] [18]. This factor sometimes appears in the 
previous membrane literature {e.g. [19]), but one gets the impression that it matters only 
in the low-stiffness regime beyond the persistence length; again, we will see that it is 
needed to get the correct correlation functions even in the more physical stiff regime. Our 
correction differs in detail from the Liouville factor because the latter is appropriate for 
conformal gauge, but in each case the motivation is the same. We will also argue that the 
usual interpretation of the Liouville factor as a Jacobian is misleading; really as argued 
above it reflects renormalization. 

Again our counterterm is important because without it, ordinary diagrammatic pertur- 
bation theory gives incorrect results due to the subtlety in the statistical measure discussed 
above. Finally we conclude in sect. 5. 

2. General Arguments 

We will consider fluid membranes, those with no internal order or structure. Real 
membranes are typically fluid at temperatures high enough to destroy in-plane order, 
but in any case such order can easily be incorporated into the argument below [3]. The 



constituents of a membrane have a preferred spacing, with a high energy cost for local 
deviations from that density. Thus, in this regime the free energy cost of a membrane 
configuration depends only on its shape. 

As mentioned in the Introduction, we can thus capture the physics of length scales 
much longer than the constituent molecule size by regarding the surface as continuous and 
using the Canham-Helfrich free energy [4] 

„.,„/d. + ^/d5(i- + i-)^ . (2.1) 

Here dS is the element of surface area, / dS is the total area of the surface, and R^^, R2 are 
the principal radii of curvature of the surface. The two coefficients /Uq, Kq are respectively 
the area cost and bending rigidity. 

The form of (2.1) is severely constrained by the requirement that the free energy 
depend only on the membrane's shape, and not on how one chooses to label points on 
the surface. While this property is manifest in (2.1), to do calculations a more explicit 
version of this equation proves necessary, in which we do choose coordinates u = (w^, u'^) 
to label points R(w) of the surface. As usual, we then get tangent vectors e^ = d^H = 
dTi/du"' at u and a metric tensor (7^^ = e^ ■ e^ at each point. The unit normal is then 
n = (e^^ X e2)/\e-^ x 62], and the curvature tensor is K^^^ = n ■ V„9^R, where V„ is the 
covariant derivative associated to the metric (7^^. Letting g = det[(7^^], [g""^] = [gab]~^^ 
and K^ = g^'^K^^^, eqn. (2.1) becomes 

n = i,, j cfu^+^ j cPu^g{K:f . (2.2) 

If we change the coordinates from -u" to u"' = u"'{u) we just relabel points in the plane. 
For an infinitesimal change u"' = u"' + e"'{u) we get a new Ti{u) = R(w) describing the 
same surface: 

R(w) = R(w) - £"(w)e„(w) . (2.3) 

The change of R is a tangential motion; conversely, tangential changes of R are not physical 
changes to the shape of the surface and ought not to be counted separately in the statistical 
sum. We thus need to supplement (2.2) by choosing a protocol for assigning a single 
coordinate system to each shape, and a procedure for discarding every surface R('u) not 
parameterized in this way. Any such choice is called a "gauge-fixing" procedure. 



A popular class of gauge- fixing procedures are the so-called "normal gauges." Suppose 
we have arranged that configurations very close to one reference surface will dominate the 
statistical sum (by stretching the membrane across a fixed frame, for instance) . We choose 
once and for all a parameterization Ro(w.) for this reference surface and compute its unit 
normal nQ{u). Then other nearby surfaces can be written as 

R('u) = Kq{u) + /i(w)no('u) (2.4) 

where the height field h at any point is uniquely defined as the normal distance from that 
point to the reference surface. We thus describe each distinct surface once when we sum 
over all functions h{u). Of course the point of this paper is that the correct definition of 
this sum is a tricky business, but let us proceed formally to obtain some general properties 
constraining the correct prescription. 

Let us consider a membrane confined to span a square, fiat frame of area A^. We 
suppose surfactant molecules can join or leave the membrane with a fixed cost in free 
energy of Hq per added area (one can easily pass to an ensemble with a fixed number of 
molecules by a Legendre transformation). The full free energy of our system Z{A^) then 
depends on the base area, and implicitly on (Uq, the stiffness Kq appearing in (2.1), the size 
a = 27r/A of the constituent molecules, and the temperature T. 

Since our frame is fiat, the equilibrium (zero-temperature) surface will be fiat too, and 
so it is convenient to work in Monge gauge, a normal gauge with Ro('U') the surface spanned 
by the frame, with Cartesian coordinates u"' running from zero to ^/A^. In addition to 
Z{A^) we will find it convenient to introduce the "effective action" r[/i]. To define T[h] 
we introduce a forcing term into 7i, 7i = 7i — / hj, and adjust j to ensure that {h{u)) is 
the desired function h.'^ Computing the partition function 

Z[A^J] = l[Dh]e-^/^e^'/^^I'^ , (2.5) 



we then let 



T[h] = -TlogZ{AsJ) + Jhj . (2.6) 



This is a formal trick to obtain correlations of height; we do not imagine j as coming from 
any pliysical force. Mathematically j = ji2d^u is a density on parameter space, so that 5/5j{u) 
is a scalar. For more details see [3]. 



The frame tension is the free energy per area of the unforced membrane, so we have 
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(2.7) 



Physically there is nothing special about Monge gauge. Even if we agree to work in 
Monge gauge there is nothing special about a reference surface Ro(w) coinciding with the 
plane of the frame; a tilted reference surface must give the same free energy. All that 
matters is the physical area of the frame. Thus the terms of T[h] involving at most first 
derivatives of h (higher derivatives are insensitive to tilting) must enter in the combination 
/ d^w ^/g, where 

9ab = ^ab + djidji 



is the induced metric and g = det (7^^ as usual. In other words, 



T[h]=T f d\{l + ^{dhf + ■■■) + 



{2.1 



where the first ellipsis involves Oilv) while the second involves more derivatives than /I's. 
We gave a more detailed argument for this conclusion in [3]. 

Whatever the precise form of the measure [Dh] appearing in (2.5), it does not depend 
on J, and so the first moment 



{h{u))^ = f[Dh]e-'^/^e^'/^^I'''h{u)=TZ 



6j{u) 



(2.9) 



Differentiating (2.6) and using (h)- = h, we find 



J{U)-TZ / e7/ N rv /^ + / eF/ N ^(^ ) =JW 
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Differentiating again and using (h) j^q = 0, we see that 



(2.10) 



lrr2 



{h{u)h{u'))^^, = Z^^T 



b'^Z 



5j{u)5j{u') 



Z-^T 



6h{u' 



5j{u) 



Sh{u'] 



5j{u) 



T 



{Zh{u')) 



52r 



5h{u)5h{u'] 



(2.11) 



Thus the inverse of the two-point function is just the quadratic part of the effective action. 
We have repeated the above well-known argument to make a point: It in no way de- 
pended upon the niceties of the measure [Dh]. Any exotic gauge-fixing or cutoff-stretching 
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factor which we may eventuaUy fold into [Dh] will not alter the above conclusion, though 
they can and will affect the relation of the quadratic part to the constant part F*-*^-*. Com- 
bining (2.11) with (2.8), (1.1) at once shows that the q^ coefficient r equals the frame 
tension r. 

Let us make a first attempt at calculating r and r. If the temperature T is small we 
can expand everything around the equilibrium configuration h — 0. The energy functional 
(2.2) then becomes [11] 

n = fj.o f d^u (1 + ^{dhf - lidhf) 

J (2.12) 

+^ / d^w [{d^hf - \{dhf{d^hf - 2{d''h)d^hd^hd^d^h] + c(/i^) . 

In this formula d^ = d/du"', d^ = 5"'' 9^ 9^ is the flat laplacian, and all indices are contracted 
with the Kronecker symbol. Thus all /i-dependence is explicitly in view in (2.12). To 
calculate Z{A^) we need a proposal for [Dh]. The simplest choice is to choose a mesh of 
grid points {u^j} on u space, spaced by a, and let 

We introduced the thermal wavelength A here to render the measure properly dimension- 
less.^ We will argue later that [Dh]^^^^^ is wrong, but for now we take it as our starting 
point. 

For low temperature we approximate Z{A^) as e~t^oAB/T ^jj^ies a gaussian integral. 
Collecting the quadratic terms of (2.12) and writing J d'^u = a^ ^ • we find 



r(o) ^ -TlogZ(A5)=^o^s + |El°g 






Now the sum is over the points q of reciprocal space. The upper limit of this sum is related 
to the constituent size a; for example, if in (2.13) we choose points on a square grid, then 
our Brillouin zone is a square of length 27r/a. Thus we recover the announced formula 

(1.2). 



In [15] and [3] this scale was taken to be the molecule size a or a constant multiple of it. 
This is legitimate when we do not care about the bare chemical potential, for example when we 
just want to dial ^o to zero physical tension. But /xo does have physical meaning, and Milner and 
Morse have argued that the factor in (2.f3) is the correct choice [f6]. 



Next we turn to the two-point function, which we provisionaUy define as 

{h{u)hm^,,^, = Z-^ j[Dh]^^^,e-^'^h{x)h{Q) . (2.14) 

This time however the constant term of (2.12) drops out, the quadratic term gives the 
leading contribution {h{q)h{—q)) q = — 5— ^ — j, and the quartic terms yield the desired 
thermal correction to {hh)Q via Wick's theorem. To shorten the formulas, in the rest of 
this paragraph we will drop all factors of A^, or set A^ = 1. Retaining only those terms 
contributing to r we get 
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Eqn. (2.15) is essentially the result of Meunier [13]. It certainly does not equal (1.2), 
notwithstanding our formal argument that r = r. 



3. Gauge Fixing 

Apparently our troubles stem from (2.13). How should we sum over all configurations? 
We will return in the next section to the question of the spacing of grid points, but even 
with this choice made correctly, eqn. (2.13) is not correct. Consider the two surfaces R, 



R+ dR in Fig. 1. According to (2.4), we describe the displacement from R to R+ dR by 
an increment 6h{u) in h{u). But this displacement is only partly physical, i.e., only partly 
normal to R, since it is along nQ(w), which is not equal to the normal n{u) at R(w). The 
actual normal displacement from R to R + dR is simply n{u) ■ nQ{u)dh{u), and we need 
to replace [Dh]^^^^^ in (2.13) by J[h][Dh]^^^^^ where 



J[h] ^ l[{n{u,^)-n,iu,^)) 



(3.1) 



y 



The Jacobian J'[h] is the first correction factor (the "Faddeev- Popov determinant") men- 
tioned in the Introduction. 
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Fig. 1: Two surfaces R, R+ dR described in the normal gauge associ- 



ated to R, 
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Our derivation of (3.1) was heuristic. Nevertheless the result agrees with a more 
formal derivation which we will give in a moment, and the heuristic discussion gives J'[h] 
a very simple geometric meaning. The formula (3.1) also agrees with the approach based 
on metrics on field space [3]. For example, taking Ro(w) to be flat we have n ■ Hq = g~^^'^, 
and JT is the Monge gauge factor found in [3] . It is not unity, contrary to ref. [14] ; while it 
does not affect the one-loop renormalization of Kq, it will be crucial for r, as shown below 
and in [3]. 

Our second derivation of (3.1) is couched in language familiar from gauge theories [20]. 
It is more precise than the motivation just given for eqn. (3.1), and again reproduces the 
answer given in [3]. We think that each derivation sheds light on (and increases our 
confidence in) the final answer. In order to form a statistical sum over distinct surfaces we 
imagine integrating over all parameterized surfaces, with a correction factor to eliminate 
overcounting. Thus, we write the partition function as 



DR] 



1 

vol[R] 

9 



-n/T 



(3.2) 



Here [-DR] just describes the possible displacement of the individual mass-points on the 
surface and vol[R] is a factor to be discussed below. All we will need to know for the 
moment is that the measure [-DR] must be coordinate-invariant, since no preferred coor- 
dinate system is given on a fluid membrane; [-DR] must also be invariant under spatial 
translations and rotations. 

To finish specifying (3.2) we need to specify the extra volume factor vol[R] appearing 
there. Starting from any surface R('u) and applying all possible reparameterizations we 
sweep out a subspace of all parameterized surfaces. The volume vol[R] of this "orbit" 
may well depend on which surface we start with, but it cannot depend on how Ii{u) is 
parameterized. 

Since it proves inconvenient to include vol[R] explicitly in (3.2), we now fix a coordi- 
nate choice by introducing two constraints on Ii{u) at every point u: 

/„(R(tx))=0 a =1,2 . 

For example, in normal gauge we require 

/„(R(tx)) ^ A-XoW-(RW-RoM) = . (3.3) 

We have introduced a factor of A~^ (see (1.3)) to make /„ dimensionless. Any length scale 
will do here since our answers will be altogether independent of our choice of /^. We now 
define the Faddeev- Popov determinant Aj.[R] via 

1 ^ Af[K]-l [DK]6[aK)] . (3.4) 

The integral is only over the space of R's which are reparameterizations of R; hence we 
also have 

vol[R] = [ [DK] . (3.5) 



To keep vol[R], and hence (3.2), dimensionless we include two factors of A~^ at every point 
in (3.5). Since every R('u) = R(0(t(,)) is related to R by a reparameterization, u"' h^ ^"(w), 
we can also write J [-DR] as J[iI>0]'D[0, R], where formally 'D[0,R] = det ^^^ ^' is a 
Jacobian factor related to the change of variables. Let us now insert (3.4) into (3.1) to get 

^ = /f^^] (^^) ^"''"'^^^^/[^] j[DQ\D[Q, R]5[/„(R(0))] . 
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Since a reparameterization Q has meaning independent of the surface R we apply it to, 
we can do the O integral last. Letting R = R(0), we use the fact that [-DR], Ti., vol[R] 
and A^ are reparameterization invariant to rewrite Z as 

Z = jiDQ] jiDK] (^] e-^[^l/^A;[R]5[/jR)]P[0,R] . 

Also, the fact that [-DR] is coordinate-invariant means that small displacements of O sweep 
out equal volumes no matter where on a given orbit they act, i.e. 'D[0,R] = ^^[O, R]. 
Renaming R as R and using (3.5) we thus get 



Z = y'[DR]e-^[^]/^A^[R]5[/„(R) 



(3.6) 



The only configurations entering (3.6) are those obeying the gauge condition, as desired, 
but now we must compute A JR] from the definition (3.4). 

While (3.6) is quite general, we will now specialize to the normal gauge (2.4) associated 
to a surface Rq. The functional delta-function in (3.6) tells us we need only to examine 
A JR] for R obeying the normal-gauge condition; the delta- function in (3.4) then tells 
us we need only study R infinitesimally close to R. Thus we write R in the form (2.3), 
R(u) = R(ti) — e"'{u)e^{u). A reasonable choice^ for [DH] near a given R is then 

[DB] ^ n(d'RK)/A') • (3-7) 

In (3.7) we imagine that whatever coordinate system we choose, a grid has been laid down 
with each point carrying constant mass. Since each constituent molecule occupies a fixed 
physical area, we can equivalently specify that each grid cell of size (Aw^, Aw^) occupy 
equal areas in 3-space. (This, in turn, typically means that the Aw" themselves are not all 
the same.) With this understanding [-DR] is indeed coordinate-invariant. 

Submitting R to an infinitesimal coordinate transformation u"' ^-^ u"' + £"(w), we see 
that at each u^j Il{u^j) moves in the tangent plane; the measure d^R(tt^j)/A^ in (3.7), 
restricted to this plane, becomes 



X-^\e,xe,\dh{u^A=X-\giu^Adhiu^A . (3.8) 



The nonlinear measure correction to be discussed in the next section will not affect our 
answer; see the Appendix. 
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We can now evaluate (3.4): 






*j 



where J^bi^) = df^{^{u + e))/de = A e^ q ■ e^ in normal gauge using (3.3). Thus 

det el ■ e^ 



A,[R]=n 



^3 



^/9 



(3.9) 



To finish our job of getting an explicit expression for the partition function, we now 
express [iI>R]5[/^(R)] appearing in (3.6) in terms of the nonredundant height variable h{u) 
appearing in (2.4). We can conveniently express an arbitrary R in terms of a height change 
plus the rest: 

R(w) =Ro(w)+/i(w)no(w)+e„o(w)i;"(w) . (3.10) 

The calculation is not quite the same as the one leading to (3.8); instead of infinitesimal 
reparameterizations e" from an arbitrary normal-gauge surface R, now we have displace- 
ments from the fixed reference surface Rq. Accordingly, instead of (3.8) we find 



d''R(w) = \ gQ{u)d^vdh . 



(3.11) 



Next, (2.4) with (3.10) gives /„(R) = A'^e^o ■ e^o^)^ or with (3.10) d^'^\f^{K)) 
A^det ~ (e^ o ' ^6 o) ' ^^'^\'^"')- We can now do the v"- integrals to get from (3.6) to 



\[\-^dh{u,;^ 



■ij 



Ae-^/^ 



where 



A = n 



u 



det (e^ ■ e^ 



Yl^-^oi^ijl 



(3.12) 



(3.13) 



«j 



It is understood that H and A are evaluated on the surface given in terms of h{u) by (2.4). 
The last step follows because det (e^ q ■ e^) = (e-,^ q x 63 g) ■ {^i x ^2) = y^Hq ■ ^/gn. Note 
that each factor in (3.12) is explicitly coordinate-invariant, in contrast with the analogous 
formula in [14]. 
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Eqns. (3.12)-(3.13) are our final form for the partition function. We see that (3.13) 
indeed does agree with our heuristic formula (3.1). Speciahzing to Monge gauge, Hq = i, 
we recover the result of [3], A = J^- g 



-1/2 



An alternate form for A can be obtained using 

eb = 9b'R = ^b,o + 9bihno)=ef,Q-he^QKf,^Q + df,hnQ , (3.14) 

where Kj^^ q is the curvature tensor of the background surface. From this we find 

^ ^ -Q det (^,,^0 -hK^b,o) ^ exp^(-la"/.a> + 0{h')) , (3.15) 



^J 



^V9 



T-J 



which again differs slightly from the formula in [14]. We see that getting the correct ^Jg 
factors in the denominator is crucial to avoid a spurious term linear in h. 

Returning to the unhappy ending of sect. 3, we now see that the Monge-gauge Faddeev- 
Popov factor 



exp 



-^E|i°g(i + (^^n 



T^^ 2 



exp 



'Y,^^m 



T ^^ 2 



(3.16) 



effectively contributes a new term to (2.12), one which is already of order T. Following 
the steps leading to (2.15), we get the corrected version 



T 



^FP ~ A^o 



d^Q 

{2-kY 



1 



^0? 



i^oQ + ^0?' 



(3.17) 



On the other hand A cannot contribute to our calculation of r since being already of order 
T we simply evaluate it at /i = 0, where it equals one. We still have the problem that 

4. Nonlinear Measure 

To see where we have erred, let us first give another derivation that r = t. Note that 
there are many different Monge gauges corresponding to projections to various reference 
planes.^ So far we have used a reference plane coinciding with the plane of the frame 



Recall that a Monge gauge is a normal gauge where the reference surface Ro in Fig. 1 is flat. 
Alternatively we could imagine two different physical surfaces related by a rotation, both viewed 
in the same Monge gauge. 
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holding the membrane, so the equihbrium configuration is h = 0, but we could instead 

choose Rq tilted by a small angle t/j. We will always work to lowest nontrivial order in t(j. In 

the new coordinate system the same equilibrium surface is now described by h'{u) = i/ju'^ 

(Fig. 2); furthermore the projected area has been foreshortened in one direction, giving 

A'^ = Aq cos?/' = Ag{l — ^V'^)- Since the two calculations compute the same quantity we 

must have that 

T[AB;h = Q] = tAb 

= T[Aj,-\i,^Aj,-h'{u)=tl,u^] (4.1) 

= TAB-T\ij''AB + r\ij^Aj^ . 

In the last step we wrote the 0{h''^) term of F in terms of r{dh'Y = rt/j'^. Thus we again 
conclude r = t, contrary to our still incomplete calculation. 




Fig. 2: Tilting the reference surface foreshortens both the base area 
and the coordinate cutoff. 

What have we done wrong? Throughout our discussion of the previous section we 
assumed that the measure [-DR], and its reduction [Dh], were invariant under changes of 
parameterization of the surface. But the choice made in (2.13) does not have this property. 
We chose the grid points {u^j} evenly spaced in the coordinates w", which has different 
meaning for different coordinate systems. In particular for our second description we have 
a total of A'^/a^ degrees of freedom, which is different from the true number A^/a^ of 
mass points! More formally, eqn. (3.7) tells us to use grid cells of constant metric area 
Au^Au'^ = a'^/y/g, while (2.13) specifies constant coordinate area Au^Au^ = a^ . We can 
easily show that this error accounts for the remaining discrepancy between r and rpp, 
eqns. (1.2) and (3.17). 

Let us denote by Fpp[A^, A^, A^; /i] the effective action calculated with the correct 
Faddeev-Popov factor J but still using the naive (coordinate-dependent) measure (2.13). 
We momentarily allow an anisotropic cutoff (A^, A^), which we write explicitly; at the end 
of the discussion we will set A^ = A^ = 27r/a. The foregoing discussion implies that we 
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may use Fpp if we remember to foreshorten the spacing of constituents (increase A^) as 
weU as the base plane (Fig. 2), so (4.1) becomes 

^ r , 27r 27r - . 
TppfAs,— ,— ,/i = 0] =rAB 

r 2.r 2.r i ^"^'^^ 

Aj, - U^Ab, (1 + W)—, — ; h'{u) = ^u' 



Tpp 



^ ^ a a 



or 



rA^ = rA, - Ir^^A^ + Ir^.^A^ + i^^^ dV^As,^^^^ ,h - 0] 



1 

Ai=A2=27r/a 

(4.3) 



since rpp is the q coefficient we calculated with the naive measure. Thus we find 

" = '■- + 91^ ■ ("> 

which really is obeyed by our calculated expressions (1.2), (3.17). 

We have thus identified the problem: the naive measure spreads mass points evenly in 
coordinate space, not in real space. For flat equilibrium surfaces {h linear in w"), the right 
prescription was easy to find: we keep the grid points uniform but change their density to 
get the RHS of (4.2). Clearly for nonflat h{u) we will have to be more clever, since here 
the cutoff is nonuniform and indeed field-dependent: the desired measure [Dh] is nonlinear 
in h. Incidentally, the same problem we face here arises again in the dynamics context [5] , 
where a naive calculation gives the two-point function controlled by our rpp. 

We hinted at the answer in the Introduction: a change in the cutoff, even a spatially- 
varying change, amounts to a local change in the number of degrees of freedom. We can 
imagine passing from the desired nonlinear [Dh] to the convenient [Dh]^^^^^ by a sort of 
decimation procedure, where we integrate e"^'^ over the discrepant degrees of freedom to 
obtain 6"^°""'^. The difference dTi = TYgg — 7i is simply a counterterm; while in general 
it is a nonlocal functional of h, still when we ask long-scale questions it may be replaced 
by a renormalizable truncation, just as we did for 7i itself. 

We therefore seek a correction factor dTi with the property that 

"[m]=e-^^['^]/^[m]„,,,," . (4.5) 

This equation is in quotes because as it stands it makes no sense: we cannot equate 
two measures with different numbers of degrees of freedom. In particular dTi is not a 
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Jacobian factor, as often asserted^, but rather is a counterterm. As such it can and will 
depend on 7i itself. Eqn. (4.5) is not a change of variables but rather an assertion about 
renormalization, namely that each side when weighted by e~^' "^ yields identical moments 
of h, at least for long wavelengths. The usual miracle of renormalization is that adjusting 
just a few coefficients in 67i suffices to make all the moments agree. 

What is 57i? We have argued that it is renormalizable. Since it arises from fluctuations 
it will be explicitly of order T, like the Faddeev-Popov correction, and so we will only need 
its form to 0{h?) in order to calculate {hh) to one loop. (In the Appendix we show how 
to do better than this.) Dimension counting says it can have at most four derivatives; 
translation invariance requires at least two. Furthermore STi should depend on h only via 
the induced metric g = 1 -\- (dh)'^, since it reflects a change of cutoff controlled by this 
metric, and in particular SH should vanish at /i = 0, since here (7^^ = 5^^ is the metric 
assumed in [-D/i] naive- Putting all these arguments together we can only have 

5n=\T^i^ I d'^u{dhf + 0{h^) , (4.6) 

where ^^ is some cutoff-dependent constant we are to flnd. (The reader who is uncomfort- 
able with the logic of this paragraph will flnd in the Appendix a more deductive derivation 
of STi obtained by matching two calculations of r[/i].) 

We remark that our STi is analogous to the "Liouville correction factor" familiar from 
the string-theory literature [17] [18]. In these papers the authors argue in analogy with the 
conformal anomaly [21] that a covariantly cut-off measure (typically the measure for the 
scale factor e*^ of an internal metric) "equals" a measure cut off by a flxed metric ((7o)a6 
times a correction: 

"[I)a]=e^^M/^[i;a]^/' , 

where 

5n[a] =T I v/^d^w [a{{daf + R^a) + ^i^e''] (4.7) 

and Rq is the gaussian curvature of g^. We are not using an internal metric, but our (4.6) 
is reminiscent of the last term of (4.7) (plus a constant). 

The pleasant feature of our approach is that flxing the coefficient /U^^ by comparison to 
a simple situation then lets us compute r[/i] in any situation. The simple situation we have 
in mind is the tilted but flat reference surface appearing on the RHS of (4.2); 57Y[/i = '^u^] 



We thank J. Distler for emphasizing this point. 
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is supposed to reproduce the effect of cutoff-stretching, which we introduced by hand via 
the last term of (4.3). Thus we choose /j^ = q ifj ^i , or 



- 1 f -^ 



log(^o?^ + '^o?' 



(aA)2 ^0^2 



+ 1- 



27rT KqQ* + ix^q 



(4.8) 



Eqn. (4.8) fixes the effective Hamiltonian to order h? . Since (4.6) is effectively just 
an additional contribution to the q^ coefficient of {hh), we see that indeed the choice (4.8) 
required by rotational invariance is also the choice which secures r = r. To calculate 
four-point correlators, or two-loop corrections to the effective couplings /Ugg, K^fj we would 
need also the 0{h^) terms of (4.6); see the Appendix. 

5. Conclusion 

If we have explained ourselves clearly the reader may find our analysis somewhat 
tautological. After all, we uncovered a paradox, only to eliminate it with a mysterious 
counterterm chosen solely for that purpose! Let us now comment on the real content of 
our analysis. 

We have studied the construction of statistical ensembles of random surfaces made out 
of fixed-size constituents. Such ensembles are good models for the long-scale behavior of 
fluid bilayer membranes, though it is straightforward to introduce in-plane order as well. 
The correct statistical weight for such an ensemble involves a subtle nonlinear measure 
[Dh] not directly suited to diagrammatic perturbation theory. We have argued that this 
measure may be replaced by a simple one if we incorporate two correction factors into an 
effective Hamiltonian. 

We have given simple geometrical interpretations to the two correction factors, as 
follows: the Faddeev-Popov term describes how a change 5h of the normal displacement 
may not itself be a normal displacement (Fig. 1). The Liouville term describes how the 
molecules of a surface, projected down to a reference plane, seem to crowd together in 
regions where the surface is tilted relative to that plane (Fig. 2). These simple pictures, 
plus some general arguments, enabled us to write down the complete form of the effective 
Hamiltonian up to terms of order T. This accuracy sufficed to show these corrections are 
crucial to secure the covariant form of the effective action argued on general grounds in 
sect. 2. Besides being satisfying, the geometrical arguments made it easier to get the right 
answers than other more formal arguments. 
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Our argument feels tautological because we used one consistency condition (r = r) to 
fix one counterterm (pi). In the appendix we sketch how to fix the rest of 67i to one loop 
accuracy, but to go farther we would need to fix its form in advance by imposing some sort 
of Ward identity, perhaps along the lines of [18]. This task we leave for future work. 
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Appendix A. Some Calculations 

In Section 4 we used rotational invariance to fix the value of part of the cutoff- 
stretching factor; here we present another approach. Our strategy is to compute the 
effective action in two different ways; demanding that the results of these calculations 
agree determines the cutoff-stretching factor to one- loop accuracy, i.e. leading nontrivial 
order in T. 

We follow the standard procedure for computing the one-loop effective action, r[h] 
[21]. We write h a-s h + dh, expand the energy functional to quadratic order in dh, and 
drop the term linear in dh. Then we do the gaussian integral over dh to get the one-loop 
contribution dr[h] to the effective action. Our key observation is that for gaussian integrals 
we don't really need a measure on field space, since in the end we just compute functional 
determinants, and for this we need only a metric on fiuctuations 6h about the given h. 
Equivalently, to one-loop accuracy we make no error if we replace the full measure [Ddh] , 
which is nonlinear in dh, by the cut off measure 



[B«i,^n^^ . (A.1) 



A 

^3,9 



where the notation "ij, (/" means that the grid points are spaced uniformly in terms of the 
distance given by the background metric g^^^ = 5„^ + djidiji, with density a~^ = (A/27r)^. 
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As we emphasize below this does not necessarily mean the points are spaced uniformly in 
any given coordinate system. The measure (A.l) is independent of dh, unlike the desired 
[DSh], but both correspond to the same metric on fluctuations about h. 

Stated succinctly, in our first calculation we cut off using g and no cutoff-stretching 
correction (though we do need the Faddeev- Popov factor as usual). 

It is a bit unsatisfying to have to use a different measure (A.l) for every background 
h. Moreover, while we succeeded in determining the free energy for various backgrounds 
h, this is not quite the same as finding the correlations of height; for the latter, the 
5/i-dependence of the true measure omitted from (A.l) will matter. So in our second 
calculation we replace [Ddh] not by [Ddh]g but by 

g-5HM/T[^^/^]^^_ (A.2) 

Here [Ddh] is defined analogously to (A.l) except that now g^ ^^ = 6^^ is the usual flat 
metric, so there is no explicit dependence on the background h: we have one measure from 
which we compute the various T\li\. Since gQ is flat, in this calculation the grid points really 
are spaced uniformly in the coordinates u"". On the other hand, now we have to expect 
a correction 671 as discussed in the text. We will find 671 to one loop by recomputing 
V\h] and comparing to the first calculation; it turns out that a single universal dTi secures 
agreement with all the different F's in the first calculation, which used a different measure 
for each h. With STi in hand we can then proceed to the calculation of {hh). 

We begin our calculations by quoting some Monge gauge formulas. We have the 
following expressions for the volume element, the inverse metric, and the mean curvature 

[11]: 

V^d^w = ^l + {dhYd'^u , (A.3) 

9 -K, l + i^dhY ' ^^-^^ 

rT^v_ -^'^ I d,d,hd^hd,h 

TTTMF + (1 + (a/i)2)3/2 • ^^-^^ 



The energy functional (2.2) becomes 



Kir 



{d^hf 2d^hd,hd,hdJ,h , {d^hd,hd,d,hf 



n = -^ d^ui :" : .„ - "',r''::r::^,r + 



v/l + (c>/i)2 (1 + (9/1)2)3/2 (1 + (9/1)2)5/2 

+ ;Uo / d^w Vl + (9/i)2 . 

19 



(A.6) 



We work to all orders in the height. 

We are now ready for the first of our two calculations. Following the above program, 
we replace h by h + 6h and retain the terms quadratic in dh. In the following we drop all 
terms that do not contribute to the frame tension renormalization, i.e. terms with more 
derivatives than /I's, though it's not much harder to retain them. We then have the 
quadratic bit of 7i, 






(A.7) 



Note that we have pulled out a factor of ^/g so that upon discretizing this integral (A.7) 
the grid points are spaced with uniform density according to the metric ^^^ (see (A.l)). 
The result of this gaussian integration is 



T„ , /KnA^a^xl 



5rga„,,,aJ/i] = -Tr^log(^ 



2Tn Jg 






(A.^ 



In this equation the symbol "Tr-" means that we evaluate the trace of the following oper- 
ator using a cutoff equivalent to the point spacing specified in (A.l). In momentum space 
this means that we find a basis of eigenmodes of the covariant laplacian A = -y=d^g'^^y/gdfjj 
then sum over all modes with eigenvalues less than (A)^. This prescription is coordinate- 
invariant; if ^^^ = d^ij then it just amounts to integrating over a circular Brillouin zone as 
usual. We have simplified (A. 8) by factoring out the constants and one power of ^~^. 

Eqn. (A. 8) is easy to work with because we can now write g"'^d^d^ as A plus terms 
with more derivatives than /I's, which we are dropping: 

^r,_,.J/.] = |tV- log ^^^^^^ - |tV- log^ + |tY- log (a^ - ^a) A-^ . (A.9) 



1/2 

which we write as exp 



To this expression we must add the contribution from the Faddeev- Popov factor fl • • - g 



_-^{\TTr-\og-g)_ 

5V\h] = -Tr, log ''Q'^'''''^^ + -Tr, log^ - -Tr, log^ + -Tr, log (^ - -^A] . 

(A.IO) 
Note that the noncovariant Faddeev-Popov factor exactly cancels with the noncovariant 
factor of g"^ from the expansion of the energy functional, leaving us with a contribution 
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to the effective action which is covariant to all orders in the height h. It is easy to evaluate 
the first trace: 

^5'°S^^^^^3^ = log^i^;^j^5:i ■ (All) 

The sum is just the total number of mass points a~^ J d'^u ^/g, which we will write as 






To make the last term of (A. 10) tractable we will now specialize to a special class of 
backgrounds h for which we know the eigenmodes of A explicitly, namely h{u) linear in 
u"". Examining this case suffices to let us extract the coefficients of all the (dh)'^ terms. ^ 
For such a ^^^ we have 

Tr^log(A2-^A)A-4 



Kq 






= M'W 7^ log iq^'''? + ^ia%r' A 



(A.13) 



After a linear transformation of integration variables we get 

^ ^ i^„ I rA , ^0,2 



= Jd\V^J^^ "^ ^°^ ^ ^' + ~^' ^ ^~' 



t2<A (27r)2 \^ Kg 

The final result for the renormalized frame tension is thus 



T [^ d^q ^ AV, 4 2^ ,A 1.N 

^ = f^o + Tr 77^-T^ log— — (koQ +^oQ ) , (A.14) 



2 7 (27r)2 ^ 27rT 

which is indeed the contribution to (1.2) due to the modes in question, and agrees with 
[15].8 

We now turn to our second derivation described above: we compute the part of the 
effective action involving the renormalized frame tension using the measure cut off by a 
fixed metric Qq „^ = 5^^. To space the grid points with uniform density according to Qq „^, 
we write (A. 7) as 

/// _ ^0 / j2„, /-— rr, -*- /a2 A^o 



n'' =-^ d'u ^Sh^= [A'-^A\Sh . (A.15) 

V9o9 V Kn 



Alternately we could simplify the sum by considering only the contribution from a thin shell 
in momentum space. 

We have corrected an error in the formula following (B.2) in [3]; it had an extra factor of ^/g 
(cf. (A.13)). 
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While the laplacian A associated to g (not Qq) again appears in (A. 15), now our gaussian 
integral is cut off using a basis of eigenmodes of the flat d^ : 



ST 



gaussian 



-Tr„ log^ = 



9o9 



^0 



Kr 



(A.16) 



To (A. 9) we must add the Faddeev- Popov factor and the cutoff-stretching factor. Both of 
these are already 0{T), so to one loop we simply evaluate them at h: 



2 so & 2nT 



9o9^ f^o ^ ^ 



(A.17) 



It is now easy to see what factor e ^'^m/'^ belongs in (A. 2). 57Y[/i] should be the 
functional of h which, when evaluated at h, reconciles (A.17) with (A. 10): 

Snlh] = - (Tr„ - Tr„ ) log '^^^^ (a^ - ^A) - -Tr„ log ( ^) . (A.18) 

^ J 2^ » 30^ ^ 2TvTg V Kq / 2 »o \9oJ 

We have written the correction factor this way to emphasize how at least part of it (the 
first term of (A.18)) is clearly a result of decimation. As in (A. 13) we can write the first 
term of (A.18) as the difference of the integral over an elliptical Brillouin zone minus the 
integral over a smaller, circular, zone with the same density of points A^ (Fig. 3). Thus 
the first term of (A.18) may be written as 



d'u 



d\ 



d^ W ^ ) log ^ {^0(9"%^^? + 1^09"'%%) (A.19) 



'BZ' (27r) J jbz 

where q is in BZ' if qaqi,g°'^ < A^ and q is in BZ if q^ < A^ 




Fig. 3: The origin of the cutoff-stretching counterterm as a difference 
between two Brillouin zones. The stretched cutoff A' equals A/cos'i/' if 
the background surface h is a, fiat plane tilted by an angle i/j relative to 
the reference plane. 
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We can recast (A. 18) as follows: 

]2^ \2^2 



- I / d^t. / ^ log ^ ('^Ol^^'^a?.)' + ^O^-'^a^) (A.20) 

-|/d^.v^/^log^ 

In this form we can readily see that the second line contributes a term to the free energy 
cancelling the contribution from Feynman diagrams, the third line cancels the contribution 
from the Faddeev-Popov determinant, and the remaining first line is covariant as required. 
From this form one can show that all the terms of STi proportional to (dh)'^ are just what 
are needed to satisfy the Ward-like identity (4.1) to all orders in the angle t/j. 
Finally, we can expand (A. 18) in terms of h; the first couple of terms are 



mh] = ^ld 



2 / d-? 

u 



{2n) 



2^2/,^ ^4 I ,, ^2^ ,, ^2 

Udhf 



\^^ A^a^(/toQ^ + ^oQ^) ^ ^ ^0? 



27rT i^oQ + fJ'oQ . 

+ 0{h^) , 

(A.21) 

in agreement with (4.8). 

To conclude our analysis, we finally compute the last correction to the q^ coefficient 
r. As mentioned earlier, the 0{h'^) term of (A.21) leads to a new correction to (hh), not 
accounted for in the steps leading to (3.17). Since this term is 0{T), it gives a new tree 
Feynman graph, which when added to the 1-loop part of (3.17) gives r = t. 
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